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ABSTRACT 



Context. The collapsar model for long gamma-ray bursts requires a rapidly rotating Wolf-Rayet star as progenitor. 

Aims. We test the idea of producing rapidly rotating Wolf-Rayet stars in massive close binaries through mass accretion and consecutive 

quasi-chemically homogeneous evolution — the latter had previously been shown to provide collapsars below a certain metallicity 

threshold. 

Methods. We use a 1-D hydrodynamic binary evolution code to simulate the evolution of a 16+15 Mq binary model with an initial 
orbital period of 5 days and SMC metallicity (Z=0.004). Internal differential rotation, rotationally induced mixing and magnetic fields 
are included in both components, as well as non-conservative mass and angular momentum transfer, and tidal spin-orbit coupling. 
Results. The considered binary system undergoes early Case B mass transfer. The mass donor becomes a helium star and dies as 
a Type Ib/c supernova. The mass gainer is spun-up, and internal magnetic fields efficiently transport accreted angular momentum 
into the stellar core. The orbital widening prevents subsequent tidal synchronization, and the mass gainer rejuvenates and evolves 
quasi-chemically homogeneously thereafter. The mass donor explodes 7 Myr before the collapse of the mass gainer. Assuming the 
binary to be broken-up by the supernova kick, the potential gamma-ray burst progenitor would become a runaway star with a space 
velocity of 27 km s -1 , traveling about 200 pc during its remaining lifetime. 

Conclusions. The binary channel presented here does not, as such, provide a new physical model for collapsar production, as the 
resulting stellar models are almost identical to quasi-chemically homogeneously evolving rapidly rotating single stars. However, it 
may provide a means for massive stars to obtain the required high rotation rates. Moreover, it suggests that a possibly large fraction 
of long gamma-ray bursts occurs in runaway stars. 

Key words. Stars: binary - Stars: rotation - Stars: evolution - Stars: mass-loss - Supernovae: general - Gamma rays: bursts 
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> 1. Introduction 



Long gamma-ray bursts are thought to be produ ced by a subset 
of dy i ng massive and possibly metal-poor stars fj akobsso n et all 
2005; lLanger & Normanll2006]: iModjaz et alj|2007l). Within the 
currently favored collapsar scenario dWooslevl [l993). the burst 
is produced by a rapidly rotating massiv e Wolf-Rayet (WR) star 
whos e core collapses into a black hole (M acFadven & Wooslevl 



stars — which locally are much rarer than O stars of solar metal- 
licity ( Langer & Norm an 2006). 



1999). While single star evolution models wi thout internal 
magne t ic fields can produce such configurations dPetrovic et alj 
2005b; Hirsch fet alJ 120051) . only models including magnetic 
fields are capable of reproducing the s l ow spins of y oung 
Galactic neutr on stars dHeger et alj|2005t fOtt et alj[2006h and 
white dwarfs (Sui js et alJ l2007l) Tdue to the magnetic core- 
envelope coupling during the giant stage. 

Various rather exotic binary evolution channels have been 
proposed to lead to lon g gamma-ray bursts dFrver et alJ 
119991: iFrver & Hegerl 120051) . supported b y the idea that long 
gamm a-ray bursts are very rare events (cf. iPodsiadlowski et all 
2004). The recent realization that long gam ma-ray bursts may 
have a bias toward s low metallicity (e.g.. iFvnbo et alj 12003; 
iFruchter et al.l 12006) may change the situation: rather than be- 
ing exotic, GRBs may simply represent massive low-metallicity 
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lYoon & Langerl d2005l) . lYoon et al.l d2006l) and 
IWooslev & Hegerl (2006) recently showed that below a 
certain metallicity threshold, very rapidly rotating single stars 
avoid the magnetic braking of the core through the so-called 
quasi-chemically homogeneous evolution: rotationally induced 
mixing processes keep the star close to chemical homogeneity, 
and thus the giant stage is avoided altogether. While these 
models are successful in producing models which fulfill all 
constraints of the collapsar model, they require very rapid initial 
rotation. The resulting number of long GRBs thus depends 
critically on the initial distribution of rotational velocities (IRF) 
of massive stars dYoon et alJl2006h . 

The question thus arises whether the quasi-chemically ho- 
mogeneous of massive stars can also be obtained in mass trans - 
ferring massive binary systems (Vanb everen & de Loord [l994). 
since in such systems t he mass gainer can b e spun-up to close 
to critical rotation (see iPetrovic et al.ll2005allb[). inde pendent of 
its initial rotation rate. While Petrovic et al. (2005b) addressed 
this question and obtained a negative res ult, they used restr ictive 
semiconvective mixing. As discussed in I Yoon et al.l (|2006), the 
semiconvective mixing efficiency is still weakly constrained, but 
most recent stellar evolution models apply efficient semiconvec- 
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tive mixing. Thus, here we readdress the q uestion, using mod els 
with efficient semiconvective mixing, as in lYoon et al.l (2006). 



2. Method 

Our stellar model is calcu lated with the same hydrodynamic stel- 
lar evolution code as in lYoon et all d2006t). This includes the 
effect of rotation on the stellar structure, transport of angular 
momentum and chemical species via magnetic torques and ro- 
tationally induced hydrodynamic instabilities. Stellar wind mass 
loss, in particular metallicity dependent Wolf-Rayet mass loss, 
and enhancem ent of mass loss d ue to rapid rotation, have been 
included as in lYoon et all d2006l) . 

The binary evolu tion physics of our code is described in 
iPetrovic et al.l (2005a b). It includes tidal coupling, mass and an- 
gular momentum transfer, and thermohaline mixing. The mass 
transfer efficiency is determined by the angular momentum bal- 
ance of the accreting star: The amount of accreted matter is lim- 
ited by the constraint that the angular momentum which it car- 
ries does not drive the ro tation of the star beyond critical rota- 
tion (Petr ovic et al.l 12005 al) . To determine the accreted angular 
momentum, the code solves the equation of motion of test par- 
ticles leaving the m ass donor into a Roche potential (cf. also 
iDessart et all2003h . 

Here, we appl y efficient semicon vective mixing; i.e., a value 
of asM = 1-0 (cf. lLanger et al.ll 1 985b is used in the calculations 
discussed below. However, the same binary model as discussed 
below was also computed with asM = 0.1 and qtsm = 0.01. 



3. Results 

We compute the evolution of a binary system with rotating and 
magnetic components of 16 Mq and 15 Mq, and an initial or- 
bital period of 5 days. We chose an early Case B system with 
an initial mass ratio close to one for two reasons. Firstly, the 
expected mass transfer efficiency for this case was about 60% 
(meaning that 60% of the transfered matter c an be reta i ned b y 
the mass gainer), ba sed on the calculation s by Wellstein (2001), 
lLanger et all d2004l) . and IPetrovic et all (l2005al) . Secondly, a 
Case B rather than Case A system was chosen to avoid synchro- 
nization after the major mass transfer phase. 

The evolution of the binary system proceeded as follows (cf. 
Table 1). The initial rotational velocity of both stars has been 
set to 230kms~ I , but both stars synchronize with the orbital ro- 
tation within about 1 Myr, to equatorial rotational velocities of 
only about 50kms _1 (cf. Figj2j. Rotationally induced mixing 
before the onset of mass transfer is thus neglig ible — in con- 
trast t o typical O stars evolvin g in isolation (Heger & Langer 
I2000t iMevnet & Maederl boOO). The initially more massive star 
ends core hydrogen burning after ~ 9.89 Myr, and Case B mass 
transfer begins shortly thereafter. It sheds about 12 M evolving 
into a ~ 4M helium star. About 1.5 Myr later, it sheds another 
~ 0.2 M Q as a helium giant, before exploding as Type Ib/c su- 
pernova. 

The mass gainer keeps about 6 M of the overflowing mat- 
ter, rendering the mass accretion efficiency to roughly 50%. 
Thereafter, it enters a phase of close-to-critical rotation, which 
induces rejuvenation and quasi-chemically homogeneous evolu- 
tion (Figs.|2]and[3]l. Its mass loss is enhanced by rotation. About 
5 Myr after the onset of accretion, the surface helium mass frac- 
tion of the mass gainer is increased to values above 60%, and 
Wolf-Rayet mass loss is assumed from then on. The star fin- 
ishes core hydrogen burning after another 3 Myr, at an age of 
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Fig. 1. Evolutionary track of the mass gainer in our 16 M + 
15 M early Case B binary model (5 d initial orbital period) in 
the HR diagram (solid line), from the zero age main sequence 
up to core carbon exhaustion. The main evolutionary phases are 
labeled by numbers (see legend). The dashed line shows the evo- 
lutionary track of a very rapidly rotating (finit/f k = 0.9) 24 M 
single star. Both stars have SMC metallicity, and undergo quasi- 
chemically homogeneous evolution (see text). 
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Fig. 2. Upper panel: Equatorial rotation velocity (dashed line) 
and critical rotation velocity (solid line) of the mass gainer of 
the computed 16 M +15 M early Case B binary sequence, as 
function of time, from the zero-age main sequence until core 
carbon exhaustion. Lower panel: Mass loss rate of the same 
stellar model, as function of time. The numbered evolutionary 
stages correspond to those given in Fig.[T]and Tab.Q] 



18.1 Myr, with a mass of 16.8 M , a surface helium mass frac- 
tion of 95 % , and rotating with -200 km s ~ 1 . 

After core hydrogen exhaustion, the mass gainer contracts 
and spins-up to critical rotation, which leads to a mass shed- 
ding of almost 2 M . During its remaining lifetime of less than 
0.5 Myr, it loses about another 2 M to a Wolf-Rayet wind. It 
ends its life as a rapidly rotating Wolf-Rayet star with a final 
mass of about 13 M , ready to form a collapsar. Assuming the 
binary broke up upon the explosion of the mass loser, the mass 
gainer would have traveled for about 7 Myr with its final orbital 
velocity of 27 km s _1 a distance of about 200 pc. 
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Table 1. Major evolutionary phases of the computed 16 M + 15 M early Case B binary sequence. The binary calculation ends 
after core carbon exhaustion of the mass loser (the primary), and the mass gainer (the secondary) is then evolved as a single star. 
We show evolutionary time, masses of both stars, mass lost from the system, orbital period, surface rotational velocities, central and 
surface helium mass fraction of the mass gainer, and orbital velocites of both stars. The abbreviations for the evolutionary phases 
are: ZAMS = zero age main sequence; ECHB = end core hydrogen burning; ICB= ignition of carbon burning; ECCB = end core 
carbon burning. The numbered evolutionary stages correspond to those given in Fig.Q]and Fig. [2] 
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Fig. 3. Evolution of the internal structure of the mass gainer of 
the computed 16 M + 15 M early Case B binary sequence, 
as function of time, from the zero-age main sequence to core 
carbon exhaustion. The time axis is logarithmic, with the time 
of core collapse as zero point. Convective layers are hatched. 
Semiconvective layers are marked by dots (red dots in the elec- 
tronic version). Gray (blue) shading indicates nuclear energy 
generation (color bar to the right of the figure). The topmost solid 
line denotes the surface of the star. 

Following [Brandt & Pods iadlowskil d 19951) we calculate the 
kick velocity necessary to unbind the binary system, under the 
hypothesis of instantaneous removal of the SN ejecta. Two ex- 
treme values correspond to the most and to the least efficient 
geometrical configuration for the supernova kick to break up 
the system. The minimum kick velocity is 52kms~ 1 and cor- 
responds to the case where the kick is aligned with the orbital 
velocity vector of the supernova progenitor. The maximum kick 
velocity necessary to unbind the system is 350 km s -1 , which is 
required if the kick is aligned to the orbital velocity vector, but 
directed backward. According to the observed velocity distribu- 
tion of radio pulsars, about 55% of pulsars have a space veloc- 
ity larger than SS Okms" 1 , while more than 98% have a velocity 
above 52 km s" 1 dArz oumania n et al 1 120021) . In order to estimate 
the chance of obtaining a runaway star out of our system we per- 
formed a Monte Carlo simulation for a randomly oriented su- 
pernova kick. According to the observe d velocity distribution of 
radio pulsars dArzoumanian et aT]|2002l) the probability for the 



Table 2. Average specific angular momentum in the CO core 
(< 7co >) for six different stellar evolution models, at the end 
of carbon core burning. The first three (labeled 'binary') cor- 
respond to the mass gainers of the computed 16 M + 15 M 
early Case B binary sequence, for three different values of the 
semiconvection parameter. The fourth corresponds to the com- 
puted 24 M Q single star with initially 90% of Keplerian velocity 
(yinit/^K = 0.9). The last two correspond to the 20 M single star 
models wi th Z-0.004 and in itially 60% and 30% of Keplerian 
rotation of YoonetaL (2006). Models in bold face are evolving 
quasi-chemically homogeneous. The specific angular momen- 
tum of the least stable orbit around a 3Mq Kerr black hole for 
these models is about 30 x 10 15 cm 2 s _1 . 
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binary system to break up by the first supernova explosion is 
about 80%. 

The same binary model as discussed above was also com- 
puted with asM =0.1 and q-sm = 0.01. The first case mentioned 
practically reproduces the results outlined above, even if the CO 
core angular momentum content of the GRB progenitor is lower 
in this case (see T ab. The second case confirmed the finding 
of iPetrovic et al.l ((2005b) that chemically homogeneous evolu- 
tion does not occur for restrictive semiconvective mixing. 

4. Comparison to single star model 

It is instructive to compare the evolution of the mass gainer of the 
binary model described above with that of a rapidly rotating sin- 
gle star of similar mass. Fig.|4]shows the Kippenhahn diagram of 
a 24 M single star with SMC metallicity with an initial rotation 
rate of 700 kms , corresponding to 90% of Keplerian ro tation 
(flinit/wK = 0.9). Modeling details are as in lYoon etall d2006l) . A 
comparison with Figf3]reveals that its evolution is almost identi- 
cal to that of the mass gainer after accretion in the binary model 
described above. This similarity is underpinned by a compari- 
son of the evolutionary tracks of both stars in the HR diagram 
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Fig. 4. Evolution of the internal structure of a single star model 
with initial mass of 24 Mq and initial rotation close to Keplerian 
("init/fK— 0.9). The evolution is shown from the zero-age main 
sequence to the end of O burning. Convective layers are hatched. 
Semi-convective layers are marked by dots (red dots in the elec- 
tronic version). The gray (blue) shading gives nuclear energy 
generation rates in log scale, as indicated on the right side. The 
topmost solid line denotes the surface of the star. 

(FiglD- Table 2 shows that also the final core angular momen- 
tum of the binary model is not significantly different from that 
of corresponding single stars. 

As a consequence, one may conclude that the binary model 
does not, from the point of view of the internal stellar evolu- 
tion, provide anything new or different from what is already ob- 
tained in rapidly rotating single stars. In particular, it can not 
be expected that the metallicity thresh old for obtaining a long 
gamma-ray burst (cf. lYoon et aTl l2006) can be significantly in- 
creased through the type of binary evolution considered here. 
While a rejuvenated accreting star is somewhat more evolved 
than a zero-age main sequence star, this difference is small and 
leads only to the avoidance of a small fraction of the mass loss 
induced spin-down during core hydrogen burning. However, it is 
to be said that single stars which rotate initially wi th 90% of their 
break -up velocity might not form in nature (cf. iMokiem et al.l 
2006, and see below). Thus, perhaps the main benefit of the mas- 
sive close binaries is just to produce very rapidly rotating O stars. 

5. Discussion 

The binary evolution model presented above shows that quasi- 
chemically homogeneous evolution may occur in mass gainers 
of low-metallicity massive early Case B binaries. The compar- 
ison of the mass gainer with a corresponding single star model 
made it clear that such binary components evolve in the same 
way as extremely rapidly rotating single stars. This confirms that 
the scenario of quasi-chemically homogeneous evolution might 
not be restricted to single stars, but may apply to the accreting 
component of massive close binaries as well. 

While we provide only one example, it seems likely that 
this scenario applies to most massive close binary compo- 
nents which accrete or gain an appreciable amount of mass; 
this ma y encompass Case A bin a ries and early Case B bi- 
naries (fPodsiadlowski et al.l 1 19921: IWellstein & Langer Il999t 
IWellstein et al.ll2.Q0lb . Case A merger are also likely contribut- 
ing to this scenario. While the merged object will have more 



mass than the initially more massive star in the binary, the prod- 
uct will be extremely rapidly rotating due to the orbital angula r 
momentum, as in the case of some blue stragglers dLivioll 19931) . 

5.1. Binaries and the distribution of rotational velocities 

The best constraint so far on the distribution of initial rotational 
velocities (IRF) comes from the recent study of young O stars 
in the SMC, mostly from the c l uster N GC 346 jMokiem et all 
I2006I) . According to lYoon et al.l d2006l) . the three most rapid ro- 
tators from the sample of 21 O stars would qualify for the quasi- 
chemically homogeneous evolution scenario, and remarkably, 
all three stars are found to be helium-enhanced. The simplest 
approach to understand those stars is to assume that they corre- 
spond to the tip of the IRF. 

However, that data of iMokiem et al.l (120061) reveals another 
interesting feature: two of the the three mentioned stars are run- 
away stars, with radial velocities deviating by 30... 70 km s~' 
from the average cluster radial velocity. While dealing with low 
number statistics, this information opens another possibility: that 
the most rapidly rotating young O stars in the SMC are products 
of binary evolution. A closer examination of the IRF derived 
by IMokiem et al.l d2006) appears to support this idea: While the 
three rapid rotators show v sin i > 290 km s : , all other stars have 

usin/ < 210kms~ 1 . 

The following hypothesis therefore seems conceivable: The 
IRF of single O stars in the SMC ends at about 210kms~ 1 
— too early to allow quasi-chemically homogeneous evolu- 
tion and collapsar formation. However, massive close binary 
evolution enhances the IRF to what we may call the appar- 
ent IRF as measured by Mokierr Tet al.l (12006b . whi ch leads to 
the re dshift dependent GRB rate as worked out by lYoon et ail 
(120061). According t o the b inary population synthesis model of 
Podsiadlo wski et all dl992l) . about 10% of all massive binaries 
might lead to a Case A merger or early Case B mass transfer, 
which is sufficient to populate the rapidly rotating part of the 
IRF of Mokiem et al. In that context, the rapidly rotating O star 
in the sample of Mokier rTet al.l (120061) which does not appear as 
runaway star could either have an undetected high proper mo- 
tion, or it could be the result of a Case A merger — where no 
runaway is produced. 

5.2. Effects from runaway GRBs 

The runaway nature of a GRB progenitor, as obtained in our ex- 
ample, has important observational consequences for both the 
positions of GRBs, and their afterglow properties. Concerning 
the afterglow, it is relevant that the medium close to a WR star 
has the density profile of a free-streaming wind, and analytical 
and numerical calculations both suggest that th e free wind of a 
single WR typically extends over many parsec dvan Marie et al.l 
2006). However, from the analysis of GRB afterglows, a con- 
stant c ircumstellar medium density has been inferred in many 
cases dChevalier & Lill2000l: iPanaitescu & Kumarl l200ll 120021: 
Chevalier et al. 20040. A possible explanation has been pro- 
posed bv lvan Marie et al] (120061) . who simulated the circumstel- 
lar medium around a moving WR star. As the GRB jet axis is 
likely perpendicular to the space velocity vector, the jet escapes 
through a region of the bow-shock where the wind termination 
shock is very close to the star. Therefore, the jet may enter a 
constant density medium quickly in this situation. 

Concerning the GRB positions, since the spin axis of the 
stars in a close binary system are likely orthogonal to the or- 
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bital plane, the observation of a GRB produced by the proposed 
binary channel is possible only if the binary orbit is seen nearly 
face on. Then the direction of motion of the runaway GRB pro- 
genitor must be orthogonal to the line of sight, allowing the pro- 
genitor, for the given space velocity, to obtain the maximum pos- 
sible apparent separation fr om its formation region. The find- 
ing of Hammer et al] d2006l) . that the nearest three long gamma- 
ray bursts may be due to runaway stars is in remarkable agree- 
ment with our scenario. While the collapsar progenitor in our 
binary model travels only 2 00 pc before i t dies, compared to 
the 400.. .800 pc deduced by Ha mmer et~aT] d2006l) . binary evo- 
lutio n resulting in higher runaway velocities are certainly possi- 
ble dPetrovic et al.ll2005a) . It remains to be analyzed whether the 
runaway scenario is compatible with the finding that long GRBs 
are more concentrated in the bright est regions of their h ost galax- 
ies than core collapse supernovae (Fru chter et alj|2006l) . 
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